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ABSTRACT 

We predict how the observed variations in galaxy populations with environ- 
ment affect the number and properties of gravitational lenses in different en- 
vironments. Two trends dominate: lensing strongly favors early-type galaxies, 
which tend to lie in dense environments, but dense environments tend to have 
a larger ratio of dwarf to giant galaxies than the field. The two effects nearly 
cancel, and the distribution of environments for lens and non-lens galaxies are 
not substantially different (lens galaxies are slightly less likely than non-lens 
galaxies to lie in groups and clusters). We predict that ~20% of lens galaxies 
are in bound groups (defined as systems with a line-of-sight velocity dispersion 
a in the range 200 < a < 500 km s"^), and another ~3% are in rich clusters 
(a > 500 km s~^). Therefore at least ~25% of lenses are likely to have environ- 
ments that significantly perturb the lensing potential. If such perturbations do 
not significantly increase the image separation, we predict that lenses in groups 
have a mean image separation that is ~ 0'.'2 smaller than that for lenses in the 
field and estimate that 20-40 lenses in groups are required to test this prediction 
with significance. The tail of the distribution of image separations is already 
illuminating. Although lensing by galactic potential wells should rarely produce 
lenses with image separations 6 ^ Q", two such lenses are seen among 49 known 
lenses, suggesting that environmental perturbations of the lensing potential can 
be significant. Further comparison of theory and data will offer a direct probe of 
the dark halos of galaxies and groups and reveal the extent to which they affect 
lensing estimates of cosmological parameters. 



1. Introduction 



Gravitational lensing of background sources by foreground galaxies offers a powerful 
probe of galaxy structure and evolution at intermediate redshifts (0.3 ^ ^ ^ 1; e.g., Keeton, 
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Kochanek & Falco 1998; Kochanek et al. 2000a). Lenses are also used to constrain the 

cosmological parameters Hq, Qm, and Q/^ (e.g., Falco, Kochanek & Munoz 1998; Helbig et 
al. 1999; Koopmans & Fassnacht 1999; Witt, Mao & Keeton 2000; and references therein). 
These applications rely on accurate models of the mass distribution responsible for the 
lensing, but the models may be complicated by contributions of the lens galaxy's environment 
to the lensing potential, especially when the lens galaxy lies in a poor group or rich cluster 
of galaxies. 

A crucial step in understanding how environment may affect lensing constraints on 
dark matter halos and cosmological parameters is determining the fraction of lens galaxies 
in groups and in clusters. The environments of most lenses are not known. Three lens 
galaxies are confirmed group members (PG 1115+080, B 1422+231, and MG 0751+2716; 
Kundic et al. 1997ab; Tonry 1998; Tonry & Kochanek 1999), and another four he in clusters 
(RXJ 0911+0551, RXJ 0921+4528, Q 0957+561, and HST 1411+5221; Young et al. 1981; 
Fischer, Schade & Barrientos 1998; Kneib, Cohen & Hjorth 2000; Munoz et al. 2000). While 
environmental effects have been included in models of some of these lenses (e.g., Schechter 
et al. 1997; Keeton & Kochanek 1997; Barkana et al. 1999; Bernstein & Fischer 1999; Chae 
1999; Keeton et al. 2000a; Lehar et al. 2000), lens environments have not been addressed 
statistically (see Keeton, Kochanek & Seljak 1997 for an initial analysis). In particular, 
previous predictions of the statistics of galaxy lenses could not examine environmental effects 
because they assumed an environment- independent galaxy luminosity function (e.g.. Turner 
1980; Turner, Ostriker & Gott 1984; Fukugita & Turner 1991; Kochanek 1993a, 1996ab; 
Maoz & Rix 1993; Wallington & Narayan 1993; Falco et al. 1998; Quast & Helbig 1999; 
Helbig et al. 1999). 

The existence of lenses in groups and clusters is not surprising. Lensing selects galaxies 
by mass and thus preferentially selects early-type galaxies, which are most common in dense 
environments (e.g., the morphology-density relation; Dressier 1980). However, recent studies 

suggest that richer environments tend to have a higher ratio of dwarf to giant galaxies 
than the field (Bromley et al. 1998b; Christlein 2000; Zabludoff & Mulchaey 2000), and 
dwarf galaxies arc much poorer lenses than giant galaxies. The morphology-density relation 
increases the chance of having a lens in a group or cluster while the dwarf-to- giant ratio 
decreases it, so the distribution of lens galaxy environments depends on how these two 
effects combine quantitatively. 

In this paper, we examine how changes in the type and luminosity distribution of galaxies 
with environment ("population variations") affect lens statistics. Using new results that 
quantify the type and environment dependence of the galaxy luminosity function (Bromley 
et al. 1998ab; Christlein 2000), we obtain the first predictions of the distribution of lens 
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galaxy environments and of statistical differences in the properties of lenses in different 
environments. We neglect contributions to the lensing potential from other galaxies and 
the extended dark halo of the group or cluster ("potential perturbations"), which introduce 
a tidal shear that distorts image configurations and a gravitational focusing that increases 
image separations (e.g., Schneider, Ehlers & Falco 1992). This assumption not only simplifies 
the analysis, but also constitutes the null hypothesis as to whether potential perturbations 
contribute (statistically) to lensing. 

The organization of the paper is as follows. In §2 we discuss the components of the sta- 
tistical calculations, including galaxy luminosity functions and the gravitational lens model. 
In §3 we use simple analytic results to highlight the trends of environmental effects in lens 
statistics. In §4 we use empirical galaxy luminosity functions to obtain quantitative predic- 
tions of the environmental effects. Finally, in §5 we offer a summary and discussion. 



We must formally describe a population of galaxies that can act as gravitational lens 
galaxies in order to compute the set of lenses they can produce. Depending on details of 

the observational sample, galaxies may be binned by type, environment, and/or redshift. 
Within each discrete bin, the distribution of galaxy luminosities is described using a smooth 
galaxy luminosity function (GLF) (f){L) such that (f){L) dL is the comoving number density 
of galaxies with luminosity between L and L + dL. The GLF is usually parameterized as a 
Schechter (1976) function with the form 



where i labels the type, environment, and/or redshift bin. The GLF is described by a 
characteristic comoving number density n*, a characteristic luminosity L*, and a faint end 
slope CKj. The differences between GLFs in different bins are characterized by differences in 
the Schechter parameters n*, L*, and ccj. 

Current surveys are not large enough to bin galaxies by type, environment, and red- 
shift simultaneously; only recently has it become possible to examine even two of the three 
quantities. Because we are interested in the effects of environment, we adopt the sample 
of galaxies from the Las Campanas Redshift Survey (LCRS, Shectman et al. 1996). This 
survey provides a large- volume, R-band selected sample of the nearby universe that offers 



2. 



Data and Methods 



2.1. Galaxy populations 
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four benefits for our study. First, it samples down to an estimated completeness limit of 
Mji = —17.5 + 5 log h, or nearly three magnitudes below M*. Second, its large volume offers 
a fair sample of environments from the field to groups and clusters. Third, its large number 
of galaxies means that galaxies can be binned by both type and environment with sufficient 
statistics to discriminate between GLFs in different bins. 

The fourth benefit is that there are two independent and complementary analyses of type 
and environmental dependences in the LCRS GLFs. Bromley et al. (1998ab, hereafter B98) 
classify galaxies using six type and two environment classes. They define six spectral types 
(which they call "clans" ) derived from a principal component analysis of important spectral 
features. The clans smoothly span the range from quiescent galaxies that have substantial 
absorption lines, 4000 A breaks, and old stellar populations (Clans 1 and 2) to star forming 
galaxies with prominent emission lines and a significant fraction of young stars (Clans 5 
and 6). B98 define two environment categories based on the local 3-d number density of 
galaxies: "high density" environments correspond roughly to groups and clusters of galaxies 
(identified using a friends-of- friends algorithm), while "low density" environments include 
only galaxies that lie outside groups. Table |l| gives the parameters for the GLFs derived by 
B98. 

The study by Christlein (2000, hereafter COO) uses a different scheme for classifying 
type and environment, placing more emphasis on environment than type. COO defines two 
galaxy types, using the equivalent width of [Oil] A3727 A to separate "emission line" (>5 
A; EL) from "non-emission line" (<5 A; NEL) galaxies. He defines environments by iden- 
tifying groups (using a friends-of-friends algorithm) and classifying galaxies by the velocity 
dispersion of the group in which they reside. Table H gives the parameters for the GLFs 
derived by COO. 

B98 and COO find similar trends in the type and environment dependence of the GLF. 
The characteristic density n* varies with both type and environment, but there is no obvious 
trend. The faint end slope a systematically steepens from early-type to late-type galaxies and 
from the field to denser environments, leading to an increase in the dwarf-to-giant ratio with 
local galaxy density. The fact that both studies reveal similar trends using different galaxy 
classification techniques suggests that the trends are robust. Zabludoff & Mulchaey (2000) 
detect a similar variation in the dwarf-to-giant ratio with environment using a completely 
different sample of galaxies. 

A drawback of the LCRS for our analysis is that the sample is drawn from the nearby 
universe {{z) ~ 0.1, Shectman et al. 1996), while most lens galaxies are found at redshifts 
between 0.3 and 1. There are, however, several independent studies of redshift evolution in 
the GLF (e.g., Lilly et al. 1995; Lin et al. 1999). We cannot use these studies directly because 
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they are too small to permit environment classification, but we can apply the inferred redshift 
trends. We adopt as our main hypothesis that the only evolution is passive luminosity 
evolution, but we also examine how number density evolution would afi^ect our conclusions. 



2.2. Lens model 

The gravitational lens model specifies what kind of lensed images can be produced by 
a galaxy of a given luminosity, type, and environment. We adopt the standard singular 
isothermal sphere (SIS) lens model because of its analytic simplicity, and because it is con- 
sistent with stellar dynamical models (e.g., Rix et al. 1997b), X-ray galaxies (e.g., Fabbiano 
1989), models of individual lens systems (e.g., Kochanek 1995; Grogin & Narayan 1996), 
and lens statistics (e.g., Maoz & Rix 1993; Kochanek 1993a, 1996a). The image separation 
produced by an SIS lens is independent of the impact parameter of the source relative to 
the lens and is given by (e.g., Schneider et al. 1992) 

where a is the velocity dispersion of the lens galaxy, and Dos and Dis are proper motion 
distances from the observer to the source and from the lens to the source, respectively. (The 
distance ratio is the same if the two distances are taken to be angular diameter distances.) 
The cross section for multiple imaging (in angular units) is 

A^\e\ (3) 

It is convenient to introduce a characteristic angular scale for lensing, 

r^87r(^)', (4) 

which is the image separation produced by an L* lens galaxy (with velocity dispersion a*) 
for a source at infinity. 

The SIS lens model omits three features that are known to affect lensing: ellipticity in 
the lens galaxy, tidal shear from objects near the lens galaxy, and gravitational focusing (or 
"convergence") from the environment. Statistically, ellipticity and shear mainly affect the 
relative numbers of two- and four-image lenses (e.g., Keeton et al. 1997), which we do not 
differentiate. Convergence from the environment can increase image separations by a few 
percent up to ~20% (e.g., Bernstein & Fischer 1999). Our use of the SIS lens model amounts 
to an assumption that environmental contributions to the lensing potential do not affect lens 
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statistics; thus our results will serve as the null hypothesis when examining whether potential 
perturbations affect lens statistics. 

To apply SIS lens models to galaxies classified by luminosity, we must convert from 
luminosity to velocity dispersion using empirical scaling relations like the Faber- Jackson 
relation for early-type galaxies and the Tully-Fisher relation for late-type galaxies. Both 
relations have the form 



L 



a 



L* \(7*J ' ^^"^ 
For early-type galaxies, combining the Faber- Jackson relation with gravitational lens statis- 
tics yields (Kochanek 1993a) 

7^4, a* = 220 ± 20 km s-\ r = 2'.'79, (6) 

which also agrees with dark matter models for the stellar dynamics of elliptical galaxies 
(Kochanek 1994). For spiral galaxies, we find that combining the R-band GLF (Tables |I] 
and ^ with the R-band Tully-Fisher relation of Sakai et al. (2000) yields 

7 = 3.6 ±0.3, a* = 100 ± 6 km s-\ r = 0'.'58. (7) 

This is an updated version of the normalization given by Fukugita & Turner (1991). The 
normalization differences between early-type and late-type galaxies are basically differences 
in mass. Early-type galaxies tend to have more mass for a given luminosity (a higher mass- 
to-light ratio) than late-type galaxies, at least within the optical radius, so they tend to have 
higher velocity dispersions and produce larger image separations. 

We assume that a* and 7 do not change with redshift. Under passive luminosity evo- 
lution, the mass scale given by cr* does not evolve. Also, high-redshift studies of the Tully- 
Fisher relation (e.g., Vogt et al. 1996; Rix et al. 1997a) and of the Fundamental Plane of 
elliptical galaxies (e.g., van Dokkum & Franx 1996; Kelson et al. 1997; van Dokkum et 
al. 1998; Kochanek et al. 2000a) suggest that 7 does not evolve. Even if the no-evolution 
assumption is not entirely justified, Mao (1991), Mao & Kochanek (1994), and Rix et al. 
(1994) have shown that lens statistics are affected by evolution only if there are dramatic 
changes in the population of early-type galaxies at redshifts z < 1. 



2.3. Statistics 



We want to compute statistical properties of both non-lens and lens galaxies. For sim- 
plicity, in this discussion we refer to the set of galaxies in a particular type and environment 
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bin i as galaxy population i. The total number density rii and luminosity density pi^i for 
galaxy population i are obtained by integrating over the luminosity function, 

POO 

rii = / ^i{L)dL = n*T{l + ai- L,^t/L*) , (8) 

roo 

PL,i = / L (PiiL) dL = n* L* r (2 + ; L,^t/L*) , (9) 

-^cut 

where the integrals include only galaxies brighter than some completeness limit Lcut, and 
they can be evaluated in terms of incomplete F functions. 

For lens galaxies, the main statistical quantity is the optical depth for a source at redshift 
Zg to be lensed by galaxy population i (e.g.. Turner et al. 1984; Fukugita & Turner 1991), 

Tiiz,) = l.JdvJdL ML, zi, Zs) ML) , (10) 

where (pi is the GLF for the galaxy population and Ai{L, zi, Zg) is the cross section for 
multiple imaging by a lens of luminosity L at redshift zi. In order to compute the number 
of lenses expected in any real survey, the optical depth would have to be modified by a 
"magnification bias" factor to account for lenses where the source is intrinsically fainter 
than the flux limit but magnified above the threshold (e.g.. Turner 1980; Turner et al. 1984). 
For the SIS lens model, magnification bias simply yields a coefficient multiplying the optical 
depth, and the coefficient is the same for all galaxy populations. Because we are mainly 
interested in the relative number of lenses produced by different galaxy populations, we can 
neglect magnification bias. 

If there is no redshift evolution in the GLF other than passive luminosity evolution, 
the optical depth and other related quantities can be computed analytically. The following 
results were originally given by Gott, Park & Lee (1989), Fukugita & Turner (1991), Kaiser 
(1991), and Kochanek (1993b); we have modified them to introduce an explicit luminosity 
cut at the completeness limit Lcut- The total optical depth is 



niz,) = T* fizs) F (1 + a, + 477! ; Le,t/A*) , (H) 

/cr*\ 4 

where r* = 16n^ rln* (^) , (12) 



c 



/Us) - ^3 L q + f],Z^yr|,)V2^^°'' ^^^^ 



H 







where Doi is the proper motion distance from the observer to the lens, th = c/Hq is the 
Hubble distance, and = 1 — — is the curvature density. Writing the optical 
depth as in eq. (0) emphasizes how the various dependences separate. The characteristic 
optical depth r* contains the characteristic number density n* and mass scale a* for the 
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galaxy population. The F function contains the shape information from the GLF (the faint 
end slope a,) and the L-a scaling relation (the Faber- Jackson or TuUy-Fisher slope 7j). 
Finally, the dimensionless factor f{zs) contains all of the dependence on cosmology and on 
the source redshift; for any fiat universe {flk = 0), /(-2s) = Dl^/ (SOrj^) (Gott et al. 1989). 
When considering a population of sources, the cosmology factor / would also include an 
integral over the source redshift distribution. 

We characterize the distribution of lenses using the differential optical depth dri/d6, 
which is equivalent (up to a normalization factor that depends on the cosmology and the 
source redshift distribution) to a probability distribution for the image separation 6. Hence 
we refer to dri/dO as the image separation distribution. It can be evaluated analytically for 
any flat cosmology, 

^ = 30r*/(^,)| (F2-2^>4 + ^^^r6) , (14) 

F„ = F[l + a,-n7ri;max(^^»/^Le,t/A*)] , (15) 

where 9 — 9/9*. Moments of this distribution can also be evaluated analytically for any fiat 
cosmology, 

_ et F(i + a. + 67-^Wi-:) 

^ ~ 2 F(l + a. + 47-;L.„,/L*) ' ^''^ 

_ 2 ie:r r (i + a. + 87r^ ; L,u.m) 

^ ~ 7 F(l + a. + 47-;W4*) ' ^ ^ 

where (■ ■ ■)i denotes an average over galaxy population i. The mean image separation is 
simply {9)i, while the standard deviation is [{9'^)i — {9)fY^'^. 

All of these results apply for a single galaxy population in a fiat cosmology. For the 
generahzation to arbitrary cosmologies, see Kochanek (1993b). For the generalization to 
multiple galaxy populations, the linearity of the optical depth means that we can simply 
sum over populations, 

Ttot = Yl'^i ^ (IS) 

i 

When computing moments of the image separation distribution, we must be careful to weight 
the populations correctly, 

(^)tot = , (20) 



9 



and similarly for {9'^)tot- 

Finally, all of these results apply to a galaxy population with no redshift evolution 
other than passive luminosity evolution. However, redshift surveys suggest that there is 
evolution in the co moving number density of late- type galaxies (Lilly et al. 1995; Lin et 
al. 1999). Adding number density evolution n*{z) to the statistical results given above is 
straightforward: evaluate t* using the local number density n*(0) and replace the cosmology 
factor f{zs) with the cosmology /evolution factor 



Note that making an explicit luminosity cut at the LCRS completeness limit Mr = 
— 17.5 + 51og/i excludes some galaxies and therefore yields an underestimate of the lensing 
optical depth. The cut is necessary because the luminosity function for any survey is unreli- 
able below the survey's completeness limit. However, the cut does not significantly affect our 
results, because the faint galaxies below it are very poor lenses. To estimate the amount of 
optical depth excluded by the cut, consider a GLF that matches an observed GLF about the 
completeness limit, but is allowed to have an arbitrary slope afaint below the completeness 
limit. Let Tbright (^faint) be the lensing optical depth due to galaxies brighter (fainter) then 
the completeness limit. Using the COO sample of NEL galaxies outside groups (see Table |^), 
we estimate Tfamt/Tbright = (0.005,0.008,0.016) when the slope below the completeness limit 
is afaint = (—0.5, —1.0, —1.5). In other words, even if the GLF is much steeper below the 
completeness limit than above the limit, the excluded optical depth is a small fraction of the 
total. 



Before making quantitative predictions about the effects of type and environment in 
lens statistics, it is instructive to use the analytic results from §2.3 to identify the general 
trends. There are two familiar effects related to galaxy type, and two new effects due to the 
environment. 

Because lensing selects galaxies by mass and early-type galaxies tend to be more massive 
than late-type galaxies, lensing has a strong type selection in favor of early-type galaxies. 
Eq. (|l5) shows that the lensing optical depth is proportional to n* (cr*)^ r(l + a + 47~^), so 
the relative number of lenses is expected to be iV(early)/A^(late) ~ 20, using the data from 
§2. Moreover, early-type galaxies tend to produce lenses with larger image separations than 
late-type galaxies. Eq. ( p!6[ ) indicates that the mean image separation (6) scales as (cr*)^. 




(21) 



3. Analytic Trends 
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so {9) for late-type lenses is only ~20% of that for early-type lenses, although this simple 
estimate may be complicated by inclination effects due to the thin disk in spiral galaxies 
(see Mailer, Flores & Primack 1997; Wang & Turner 1997; Keeton & Kochanek 1998). Both 
of these effects are known from previous calculations of lens statistics (e.g.. Turner et al. 
1984; Fukugita & Turner 1991). They are also consistent with the data, as most of the more 
than 50 known lenses are produced by early- type galaxies (e.g., Kochanek et al. 2000a), 
with only four likely cases of lensing by a spiral galaxy (B 0218+357, Browne et al. 1993; B 
1600+434, Jaunsen & Hjorth 1997; PKS 1830-211, Wiklind & Combes 1996; 2237+0305, 
Huchra et al. 1985). Also, the smallest known image separation is produced by a face-on 
spiral galaxy {6 = 0'.'33 for B 0218+357). The type dependence of the number density of 
lenses and the mean image separation leads to an environmental dependence because of the 
morphology-density relation, as we show in the next section. 

The two new environmental effects arise from the environmental dependence of the faint 
end slope a of the GLF. First, because the faint end slope appears to systematically change 
with environment (see Tables |I] and ||, and Zabludoff & Mulchaey 2000), the distribution 
of image separations should vary with environment. Figure ^ shows the image separation 
distribution dr/dO for different values of a. As a decreases, the distribution becomes more 
skewed toward small separations, and {9) decreases. Physically, when a is smaller the galaxy 
population has a higher ratio of dwarf galaxies to giant galaxies, which leads to a higher 
fraction of small-separation lenses. 

Second, changes in a with environment affect the numbers of lens and non-lens galaxies 
differently. Consider Figure^, which shows schematic representations of groups with a = 0.5 
and —1.5, normalized to have the same number of galaxies down to M* + 3. If we pick a 
galaxy at random, it is equally likely to come from either group. However, if we pick a lens 
galaxy at random, it is far more likely to come from the a = 0.5 group than the a = —1.5 
group, because the latter group is dominated by dwarf galaxies. Considering other toy 
models, such as groups normalized to have the same lensing optical depth, leads to a general 
conclusion: other things being equal (i.e., in the absence of a morphology-density relation), 
lens galaxies are less likely than non-lens galaxies to be found in environments with a large 
dwarf-to-giant ratio. The effect is probably weaker in observed galaxy populations than in 
our toy models with their rather extreme values of a. 

4. Quantitative Results 

We now evaluate lens statistics using empirical GLFs to quantify the effects of the 
morphology-density relation and the dwarf-to-giant ratio on lens statistics. In §4.1 we present 
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results using the full set of type and environment bins defined by B98 and COO, while in 
§4.2 we use coarser binning with only two types (early- and late-type galaxies) and three 
environments (the field, poor groups, and rich clusters). In §4.3 we compare our results to 
the data. In §§4.1-4.3 we assume no evolution in the comoving number density of galaxies 
with redshift, but in §4.4 we discuss how evolution would affect our results. Finally, in §4.5 
we discuss the effects of possible incompletenesses in the Las Campanas Redshift Survey. 



4.1. Fine binning 

Figures ^ and ^ show the image separation distributions dTi/dO for the type and environ- 
ment bins defined by B98 and COO, respectively. The distributions are computed assuming 
a source redshift Zg = 1 and a cosmology with VLm = 0.3 and VLiy = 0.7. If we were to change 
the cosmology or the source redshift, or even to allow a distribution of source redshifts, the 
only effect would be to change dTi/dO by a multiplicative factor that is identical for all type 
and environment bins (see eq. |Tl|). The relative contributions of the different bins to the 
lensing optical depth are independent of the cosmology and source redshift distribution. 

The results illustrate the trends identified in §3. First, for galaxies in a particular envi- 
ronment, both the amplitude and the peak location of the distributions are larger for early- 
type galaxies than for late-type galaxies; relative to early-type galaxies, late-type galaxies 
have a much smaller lensing optical depth and produce smaller image separations. Second, 
for galaxies of a particular type, the lens distributions in denser environments are more 
dominated by small image separations. This effect occurs because denser environments tend 
to have a steeper GLF and hence a higher dwarf-to-giant ratio (B98; COO; Zabludoff & 
Mulchaey 2000), which means a higher fraction of galaxies that produce small separation 
lenses. These trends are seen in both the B98 and COO samples, which suggests that they 
are not overly sensitive to details of the type and environment classification schemes. 



4.2. Coarse binning 

The type classification of B98 and the environment classification of COO both yield finer 
differentiation than we seek for characterizing the type and environment distributions of lens 
galaxies. We confine the B98 sample to two type bins by defining an early-type sample using 
Clans 1-3 and a late- type sample with Clans 4-6. We reclassify the COO sample into three 
environment bins: the field, poor groups, and rich clusters. Studies of groups suggest that 
many systems with velocity dispersions a ^ 200 km s~^ are either spurious or similar to 
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the Local Group, which is bound but not yet viriahzed (Diaferio et al. 1993; Zaritsky 1994; 
Zabludoff & Mulchaey 1998). Hence, to be conservative, we define a "field" bin that contains 
both galaxies not in groups and galaxies in systems with a < 200 km s~^, a "group" bin that 
contains galaxies in systems with 200 < a < 500 km s~^, and a "cluster" bin that contains 
galaxies in systems with a > 500 km s"^. As we demonstrate below, our main conclusions 
are robust to the effects of changing the definitions of the coarse bins. Note that we compute 
statistical quantities for the full set of type and environment bins, and then combine the 
results into the coarse bins using eqs. (|l8D -(^); we do not try to define a GLF for each 
coarse bin. 

Figure |^ shows the relative contributions of the coarse type and environment bins to 
the sets of non-lens and lens galaxies. We emphasize that the results for non-lens galaxies 
are empirical results from observed galaxy samples, while the results for lens galaxies are 
predictions based on the SIS lens model. The most dramatic result is the familiar type 
difference between non-lens and lens galaxies. Late-type or EL galaxies account for 66% and 
68% of non-lens galaxies in the B98 and COO samples (respectively), but only 4% of lens 
galaxies in either sample.[] Again, these results reiterate the well-known fact that lensing 
selects galaxies by mass and overwhelmingly favors early-type galaxies. 

The environment fractions are given in Table ^. The morphology-density relation is 
evident in the environment fractions for non-lens galaxies; it can also be seen in the NEL 
fraction, which for the COO sample is (0.30, 0.39, 0.52) in (the field, groups, clusters). Acting 
alone, the morphology-density relation would imply that lens galaxies should be more likely 
than non-lens galaxies to lie in dense environments. But lensing is also affected by changes 
in the dwarf-to-giant ratio with environment. It turns out that the two effects nearly cancel, 
and the environment fractions are very similar for non-lens and lens galaxies. In fact, the 
changing dwarf-to-giant ratio is the slightly stronger effect. Hence lens galaxies are slightly 
less likely than non-lens galaxies to lie in dense environments, although the differences are 
small and may be hard to detect. The similarity between the results computed with the B98 
and COO samples again suggests that systematic differences in classifying galaxy types and 
environments have little effect on the relative contributions of different environments to lens 
statistics. 



^Changing the type classification of the B98 Clan 3 and 4 galaxies would of course change the type 
fractions. Reclassifying Clan 3 galaxies as late-type would increase the fraction of late-type non-lens galaxies 
to 85% and the fraction of late-type lens galaxies to 11%. Conversely, reclassifying Clan 4 galaxies as early- 
type would reduce the non-lens late-type fraction to 43% and the lens late-type fraction to 1%. The effects 
are smaller for lens galaxies than for non-lens galaxies because late-type galaxies contribute a small fraction 
of lenses to begin with. 
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Quantitatively, we predict that 22-27% of lens galaxies are in high density, bound 
systems like groups and clusters — and that the majority of these are in poor groups rather 
than rich clusters. These results are consistent with known lenses, where the sample of more 
than 50 lenses includes four lens galaxies in clusters and three in confirmed groups, but most 
lens galaxy environments remain undetermined. Our predictions yield a lower limit on the 
fraction of lens systems where the lensing potential is probably perturbed by elements other 
than the lens galaxy. If the lens galaxy lies in a bound group or cluster, it is very likely that 
other group galaxies and an extended dark matter halo contribute shear and/or convergence 
to the lensing potential. Even if the lens is not in a bound system, however, there may 
still be neighboring galaxies or structures along the line of sight that noticeably perturb the 
lensing potential. 

The other trend identified in §3 is that the mean image separation depends on the faint 
end slope a, which varies with environment. Table ^ gives the predicted mean and standard 
deviation of the image separations for lenses in different environments. We predict that the 
mean separation for lenses in groups is smaller than the mean for lenses in the field, under 
the assumption that the environment does not significantly affect the image separation (see 
§2.2). The difference in the means is considerably smaller than the standard deviation, so it 
would take 20-40 group lenses to detect the difference even at the la level (assuming that 
for N items with standard deviation a the uncertainty in the mean is a/VW). Given good 
statistics, however, comparing the mean separations for group lenses and field lenses would 
provide an excellent test of environmental contributions to lensing. Specifically, an extended 
dark halo in a group containing a lens can provide extra gravitational focusing that increases 
the image separation beyond that produced by the lens galaxy alone (e.g., Falco et al. 1985). 
Thus if observed lenses reveal (^) (groups) > (^) (field) contrary to our prediction, it would 
provide direct evidence that lenses in groups are significantly affected by group dark halos. 



4.3. Comparison with data 

Figure § shows our predictions for the total image separation distribution compared 
with the data for 49 known lenses from Kochanek et al. (2000b). Such comparisons are 
sensitive to the source redshift distribution and to cosmological parameters (e.g., Kochanek 
1996ab; Falco et al. 1998; Helbig et al. 1999; and references therein). With the SIS lens 
model, these dependences appear only in a normalization factor that does not affect the 
shape of the separation distribution (see eq. [TT]), so we avoid complications by normalizing 
the theoretical curves to 49 total lenses. 

Ideally, comparisons between theory and data include a careful account of selection 
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biases, but that is impossible here because Figure ^ includes lenses from many different 
surveys as well as serendipitous discoveries. Thus it is somewhat surprising to find that the 
theoretical curves agree quite well with the data. A K-S test (e.g., Press et al. 1992) cannot 
distinguish between the data and any of the models. The agreement has little to do with 
the changes in galaxy populations with environment. Using the environment-independent 
LCRS GLFs given by Lin et al. (1986) yields a separation distribution very similar to the 
curves shown in Figure ^ The agreement is also insensitive to environmental perturbations 
of the lensing potential, which our models neglect. The largest known convergence from 
a group or cluster enclosing a lens increases the image separation by ~20% (see Bernstein 
& Fischer 1999; Romanowsky & Kochanek 1999). Randomly choosing 25% of the lenses 
in Figure ^ and increasing their image separations by 20% produces a distribution that is 
statistically indistinguishable from the original distribution. In other words, examining the 
image separation distribution without knowing the environments of the lenses is a poor way to 
test whether environments affect the lensing potential (except in the tail of the distribution, 
as explained below). A better test is to determine the environments and then compare lenses 
in groups with lenses in the field, as discussed in §4.2. 

Nevertheless, Figure ^ does offer three conclusions. First, regardless of the selection ef- 
fects that are at work, it appears that no particular class of lenses is substantially missing. For 
example, even though finite resolution might bias surveys against finding small-separation 
lenses, it appears that these lenses are not significantly underrepresented. Second, Figure ^ 
provides a reassuring consistency check: our addition of population effects to lens statistics 
has not degraded the agreement between theory and data. It has not substantially improved 
the agreement, either, because the image separation distribution is not very sensitive to 
environmental effects for the reasons just cited. 

The third point is that there are two lenses with 9 > 6", RXJ 0921+4528 and Q 
0957-1-561, that lie in the tail of the image separation distribution. The COO model shown in 
Figure |^ predicts only 0.04 lenses with 6 > 6" and 1.39 lenses with 6 > 4". Part of the expla- 
nation for these two unlikely lenses is that we have neglected environmental contributions 
to the lensing potential. In Q 0957-1-561 the lens galaxy is the brightest galaxy in a cr ~ 700 
km s~^ cluster, and the convergence from the cluster is thought to contribute ~20% of the 
image separation (e.g., Bernstein & Fischer 1999; Romanowsky & Kochanek 1999). In RXJ 
0921+4528, the lens galaxy appears to be in an X-ray cluster that may likewise contribute 
to the large separation (Munoz et al. 2000). These lenses suggest another test of our null 
hypothesis that potential perturbations do not affect lens statistics. They already indicate 
that convergence from the environment can be important when the lens lies in a cluster. As 
the data improve, it will be possible to extend this result to lower mass environments like 
groups and determine the importance of including environmental contributions to the lens 
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model in applications of lens statistics. 

4.4. Evolution effects 

Our results so far have been obtained under the assumption that the comoving number 
density of galaxies does not change with redshift. However, redshift surveys (e.g., Lilly 
et al. 1995; Lin et al. 1999) suggest that there is number density evolution, at least in the 
population of late-type galaxies. (The surveys also imply evolution in the luminosity density, 
but lensing is insensitive to the component of luminosity density evolution that is due to 
passive luminosity evolution; see §2.) To quantify number density evolution, Lin et al. (1999) 
introduce the parameter P defined by 

n{z) = n(0) IQ""-^^'. (22) 

In the redshift range 0.12 < ^ < 0.55, Lin et al. (1999) find P ^ 3 ± 1 (la) for late-type 
galaxies and essentially no number evolution in early-type galaxies. 

We can include the effects of number evolution in the lens statistics by replacing the 
cosmology factor f{zs) in eq. (|lT]) with the cosmology/evolution factor /evoi(^s) from eq. (|2T|), 
which is equivalent to multiplying the optical depth by an evolution factor femiizg)/ f{zs)- 
Table ^ shows that this factor increases strongly with both the number evolution parameter 
P and the source redshift Zg. With higher P the evolution is more rapid, and with higher 
Zs the increase in number density between the observer and the source is larger. Strong 
number evolution (P ~ 3) can increase the optical depth for late-type galaxies to lens a 
distant source [zg ~ 3) by more than an order of magnitude. The effect is unlikely to be this 
strong, however, for two reasons. First, fewer than 10% of lenses have sources beyond Zg = 3. 
Second, the strong evolution found by Lin et al. (1999) was measured in the redshift range 
0.12 < 2; < 0.55; there are no current grounds for extrapolation to 2; ~ 1 or beyond, especially 
using an exponential function that may over-estimate the evolution at high redshifts. 

Consider a rather dramatic number evolution factor of fevo\{zs) / f (zg) = 10 for late-type 
galaxies and no number evolution for early-type galaxies. Using the COO luminosity functions, 
we then find that late-type galaxies account for 32% of lens galaxies — the rapid number 
evolution in late-type galaxies simply increases the number of late-type lens galaxies. The 
mean image separations correspondingly decrease: (9) = (1'.'33, 1'.'29, 1'.'03) for (field, group, 
cluster) lenses. Nevertheless, the environment fractions are largely unchanged: 79% of lens 
galaxies in the field, 17% in groups, and 3% in clusters. (Compare with Tables |^ and ^ for the 
results without number evolution.) The environment distributions for early- and late-type 
lens galaxies are simply not very different, so changing the relative fraction of lens galaxy 
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types does not significantly change the net distribution of lens galaxy environments. 

4.5. Incompleteness effects 

All galaxy redshift surveys are incomplete with either explicit or implicit magnitude 
and surface brightness limits. Lin et al. (1996) estimate a completeness limit of Mr = 
— 17.5 + 51og/i for the LCRS, and we have restricted our analysis to galaxies brighter than 
this limit. We now ask whether any remaining incompletenesses brighter than this limit 
could affect the lens statistics. 

From a direct comparison of the R-band LCRS with the B-band CfA2 redshift survey, 
Huchra (1999) argues that the LCRS underestimates the number of galaxies with Mr <^ 
— 16 + Slog/i by a factor of ~4 and that most of the "missing" galaxies are low surface 
brightness, late-type, emission line galaxies. Although most of the suggested incompleteness 
occurs at magnitudes fainter than we consider (and fainter than the LCRS completeness 
limit) and is derived from a perilous comparison of surveys in different bandpasses, we use 
it to estimate the effects of incompleteness on our analysis. 

In §2.3 we argue that even if the GLF below the completeness limit is much steeper than 
above the limit, the optical depth from galaxies below the limit is a tiny fraction of the total. 
Hence the behavior of the GLF below the completeness limit has essentially no effect on our 
results. To check our results further, we consider the effects if the LCRS underestimates the 
number of galaxies even at magnitudes brighter than the completeness limit. We increase the 
number of galaxies in the Am = 1 mag bin above the LCRS completeness limit by the factor 
of ~4 estimated by Huchra (1999). Our approach is conservative, because we apply the 
changes at magnitudes significantly brighter than where Huchra (1999) proposes most of the 
incompleteness. If the incompleteness applies to late-type galaxies in all environments, we 
find that it changes the lens environment fractions by <0.001 compared with Table ^, and it 
decreases the mean image separations by ~0'.'02 compared with Table |^. If the incompleteness 
applies only to late-type galaxies in the field (because of the morphology-density relation), 
it changes the lens environment fractions by ~0.002 and the mean image separation in the 
field by ~0'.'02. In other words, because it mainly applies to faint galaxies that are poor 
lenses, incompleteness has a negligible effect on our results. 
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5. Conclusions 

We have studied how changes in galaxy populations with environment affect gravita- 
tional lens statistics, quantifying the results with data from the Las Campanas Redshift 
Survey (Shectman et al. 1996). Two effects determine how the distribution of lens galaxy 
environments differs from the distribution of normal galaxy environments. First, the like- 
lihood that lens galaxies lie in groups and clusters is enhanced by the morphology-density 
relation (e.g.. Dressier 1980) and the fact that lensing tends to select massive early-type 
galaxies. Second, it is diminished by the systematic increase in the fraction of dwarf galaxies 
relative to giant galaxies in dense environments (e.g., Bromley et al. 1998b; Zabludoff & 
Mulchaey 2000; Christlein 2000). The two effects nearly cancel, and the fraction of gravi- 
tational lens galaxies in a particular environment is very similar to the fraction of non-lens 
galaxies found in that environment. As a second order effect, lens galaxies are slightly less 
likely than non-lens galaxies to be found in poor groups and rich clusters. Quantitatively, we 
expect ~20% of lens galaxies to be found in poor groups (defined as systems with velocity 
dispersions in the range 200 < a < 500 km s~^), and another ~3% of lens galaxies in rich 
clusters (defined as systems with a > 500 km s~^) . 

Thus we predict that for ;^25% of lenses the lensing potential may include significant 
contributions from objects other than the lens galaxy. This result is only a lower limit for 
three reasons. First, bound groups may exist below the velocity dispersion threshold of 200 
km s~^ that we have imposed. Second, a lens that does not lie in a group or cluster can still 
be perturbed by neighboring galaxies. Finally, lensed images may be perturbed by large- 
scale structure along the line of sight. External perturbations have been included in models 
of some individual lenses (e.g., Grogin & Narayan 1996; Schechter et al. 1997; Lehar et al. 
2000), but they need to be added to applications of lens statistics such as limits on VLm and 
Qa (e.g., Kochanek 1996a; Falco et al. 1998; Helbig et al. 1999). To do so, the statistical 
analysis of environmental shear by Keeton et al. (1997) must be updated with knowledge 
of the spatial distribution of galaxies and dark matter in poor groups and rich clusters, not 
to mention correlations among member galaxies and between galaxies and the diffuse dark 
halos in groups and clusters. The improved analysis must also include lensing effects that 
arise from the non-linear interaction between ellipticity in the lens galaxy and shear from 
the environment (see Keeton, Mao & Witt 2000b). 

Our calculation provides the starting point for the improved analysis in two ways. First, 
we have assumed that environmental perturbations of the lensing potential do not affect lens 
statistics. Subsequent analyses can take our results as the null hypothesis to test whether 
potential perturbations are important. Second, we have presented empirical tests to deter- 
mine whether potential perturbations are important statistically. If they are not important. 
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we predict that the mean image separation for lenses in groups is smaller than the mean 
separation for lenses in the field, although it may take 20-40 group lenses to test this pre- 
diction with significance. A contrary empirical result would indicate the presence of extra 
gravitational focusing, not included in our models, from matter in groups and clusters that 
is not associated with the lens galaxy. A more immediate version of this test comes from the 
existence of two large-separation lenses {9 > 6", RXJ 0921-1-4528 and Q 0957-1-561), which 
is far more than predicted by our models. In fact, both lenses appear to lie in clusters, 
and in Q 0957-1-561 the cluster is thought to contribute ~20% of the image separation (see 
Bernstein & Fischer 1999; Romanowsky & Kochanek 1999). 

Out of the current sample of about 50 lenses, we expect about 10 lenses in groups and 
about two lenses in clusters. At present, the environments of most lenses have not been 
determined. Four lens galaxies appear to lie in clusters (RXJ 0911-1-0551, RXJ 0921-1-4528, 
Q 0957-1-561, and HST 1411-1-5221), and another three in spectroscopically confirmed groups 
(PG 1115-f080, B 1422+231, and MG 0751+2716). Once more lensing groups are found, 
they will provide a sample of groups at redshifts 0.3 ^ ^ ^ 1, which can be compared 
with existing samples of nearby groups (e.g., Zabludoff & Mulchaey 1998) and of distant 
clusters (e.g., Stanford, Eisenhardt & Dickinson 1998) to study galaxy evolution in different 
environments. The existence of lenses in the groups will also permit direct studies of the 
diffuse dark matter in poor groups of galaxies. 
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Table 1. B98 Luminosity Function Parameters 



Clan 


Environment 


a 


M* -51og/i 


n* 










(10-3 /i3 MpC-3) 


1 


low density 


1.10 


-20.06 


0.18 




high density 


0.20 


-20.48 


0.21 


2 


low density 


0.17 


-20.10 


5.36 




high density 


-0.39 


-20.39 


2.61 


3 


low density 


-0.09 


-19.81 


7.79 




high density 


-0.58 


-20.00 


2.97 


4 


low density 


-0.65 


-19.88 


7.01 




high density 


-0.61 


-19.78 


3.10 


5 


low density 


-1.05 


-19.80 


3.70 




high density 


-1.61 


-20.39 


0.49 


6 


low density 


-1.94 


-20.09 


1.38 




high density 


-1.93 


-20.14 


0.56 



Note. — The luminosity function parameters derived by B98 for 
their six spectral type "clans" and two environments. 
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Table 2. COO Luminosity Function Parameters 



Galaxy 


Group a 


a 


M* - 51og/i 


n* 


Type 


(km s"-^) 






(10-^ /i^ Mpc-3) 


NEL 


other 


-0.079 


-20.23 


49.31 




50 


0.202 


-20.04 


2.58 




150 


0.070 


-20.16 


4.92 




250 


-0.345 


-20.40 


5.68 




350 


-0.334 


-20.34 


6.54 




450 


-0.596 


-20.59 


2.63 




750 


-0.909 


-20.84 


2.54 


EL 


other 


-0.908 


-20.14 


63.74 




50 


-0.870 


-20.18 


3.23 




150 


-0.958 


-20.30 


4.84 




250 


-0.932 


-20.18 


6.09 




350 


-1.304 


-20.49 


3.39 




450 


-0.708 


-20.03 


3.02 




750 


-1.025 


-20.17 


2.58 



Note. — The luminosity function parameters derived by COO 
for two galaxy type and seven environment bins. One environment 
bin ("other") contains all galaxies that do not lie in groups. The 
remaining bins classify galaxies by the velocity dispersion a of the 
group in which they reside. The five environment bins labeled 50, 
150, 250, 350, and 450 km s"^ are 100 km s"^ wide and labeled by 
the central value. The bin labeled 750 km s~^ includes the range 
500 < (J < 1000 km s~^. Typical uncertainties are ~0.2 in a and 
~0.2 mag in M*; see COO for details. 
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Table 3. Environment Fractions 





B98 


Sample 




COO Sample 




Galaxy Type 


Low Density 


High Density 


Field 


Groups Clusters 


Early-type / Non-lens 


0.636 


0.364 


0.700 


0.226 


0.074 


Lens 


0.728 


0.272 


0.781 


0.186 


0.033 


Late-type / Non-lens 


0.706 


0.294 


0.796 


0.171 


0.033 


Lens 


0.733 


0.267 


0.821 


0.149 


0.030 


All / Non-lens 


0.682 


0.318 


0.765 


0.189 


0.046 


Lens 


0.728 


0.272 


0.783 


0.185 


0.033 



Note. — The fraction of non-lens and lens galaxies of different types that are ex- 
pected in each environment, computed for the samples defined in the text. We use 
Monte Carlo simulations to estimate the statistical uncertainties in the environment 
fractions due to measurement uncertainties in a and M* (given by COO) and a* (see 
§2.2). We estimate that the uncertainties are < (0.040, 0.037, 0.023) for non-lens galax- 
ies in (the field, groups, clusters), while they are <0.008 for lens galaxies in all envi- 
ronments. 

These results are robust to changes in the definitions of the type and environment 
bins. The B98 fractions are essentially unchanged if we change the type classification 
of Clans 3 and/or 4. For the COO sample, if we increase the lower velocity dispersion 
threshold for groups to a = 300 km s^^, the (field, group, cluster) fractions become 
(0.839,0.115,0.046) for non-lens galaxies and (0.853,0.114,0.034) for lens galaxies. 
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Table 4. Mean Image Separations 



Sample 


Environment 


(0) n 


n 


B98 


Low Density 


1.83 


1.04 




High Density 


1.57 


0.95 


COO 


Field 


1.76 


1.02 




Groups 


1.61 


0.96 




Clusters 


1.30 


0.85 



Note. — The predicted mean {9) and standard devia- 
tion ag of the image separation distributions for lenses in 
different environments. Monte Carlo simulations indicate 
that the statistical uncertainties in (6) are ~0'.'3, due almost 
entirely to uncertainties in a* (see §2.2). However we are 
more interested in the difference (^) (field) — (6') (groups), 
where the uncertainty is only 0''09 and is due primarily 
to uncertainties in a. (In the difference quantity, a* fac- 
tors out into a constant scale factor.) Larger surveys such 
as the Sloan Digital Sky Survey (e.g., Gunn & Weinberg 
1995) should reduce this uncertainty to ;^0'.'04 by placing 
better constraints on a. 
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Table 5. Number Evolution Factors 





Zs = 1.5 


2.0 


2.5 


3.0 


P = 1.5 


2.51 


3.28 


4.23 


5.44 


2.0 


3.53 


5.17 


7.52 


11.00 


2.5 


5.05 


8.41 


14.06 


23.90 


3.0 


7.35 


14.10 


27.56 


55.57 



Note. — The number evolution factor fcm\{zs) / f{zs) 
computed for different values of the source redshift Zs and 
the Lin et al. (1999) number evolution parameter P, as- 
suming VLm = 0.3 and f^A = 0.7. 
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Fig. 1. — Image separation distributions dr/dO for environments with different values of the 
faint end slope a of the galaxy luminosity function. All of the curves use a Faber- Jackson or 
TuUy-Fisher slope 7 = 4 and are normalized to have unit area. The mean image separation 
{9) for each curve is indicated. 
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a = 0.5 a = -1.5 




M - M* M - M* 



Fig. 2. — Schematic representations of groups with a — 0.5 (left) and a — —1.5 (right), 
normahzed to have the same number of galaxies down to M* + 3. The luminosity functions 
are shown at the bottom. For each group, galaxies are chosen randomly from the luminosity 
function, placed randomly inside the dotted circle, and drawn with a circle whose area is 
proportional to the galaxy's luminosity and lensing optical depth and whose diameter is 
proportional to the image separation. (The random placement of galaxies is not realistic, 
but it is done for illustration only and does not affect our calculations.) The points on the 
luminosity function indicate the galaxies. 
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Fig. 3. — Image separation distributions for the six clans and two environments defined by 
B98. Results are shown for a source at redshift Zs = 2 in a. cosmology with Qm = 0.3 and 
= 0.7. Each panel contains a single clan; the sohd and dotted curves denote low and 
high density environments, respectively. Note that the vertical scales differ between the two 
rows. The vertical dashed lines indicate the peaks of the distributions. The two values of T4 
in each panel give the optical depth (in units of 10~^) for low and high density environments, 
respectively. We treat Clans 1-3 as early-type galaxies and Clans 4-6 as late-type galaxies; 
if we were to reclassify Clan 3 as late-type galaxies, the results would look very similar to 
the Clan 4 panel. 
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Fig. 4. — Image separation distributions for the two types and seven environments defined 
by COO, again for Zs = 2, Qm = 0.3, and = 0.7. Each panel contains a single environ- 
ment; the solid curves denote NEL galaxies and the dotted curves EL galaxies. The vertical 
dashed lines indicate the peak of each distribution. In each panel, the two values of give 
the optical depth (in units of 10~^) for NEL and EL galaxies, respectively. Note that the 
vertical axis scale for the "other" panel differs from that in the group panels. For the lu- 
minosity/mass conversions, we treat NEL galaxies as early-type galaxies and EL galaxies as 
late- type galaxies. 



B98 non-lens 



COO non-lens 




Fig. 5. — Pie charts giving the relative contributions of the different type and environment 
bins to the set of non-lens galaxies (top) and the set of lens galaxies (bottom). Results 
are shown for the B98 (left) and COO (right) samples defined in the text. For the COO lens 
sample, the EL/clusters bin has a fraction 0.001 and is too small to be seen. 
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Fig. 6. — The histogram shows the distribution of image separations for 49 known lenses, 
in bins of width 0'/2. The sohd curve shows the net distribution for the COO sample (a sum 
of the curves shown in Figure 4). The three dashed curves show the net distributions for 
three different models based on the B98 sample. The middle curve represents the model 
discussed in the text, in which Clans 1-3 are treated as early- type galaxies and Clans 4-6 
as late- type galaxies (a sum of the curves in Figure 3). In the higher-peaked curve Clans 
1-4 arc treated as early-type galaxies, while in the lower-peaked curve only Clans 1-2 are 
treated as early-type galaxies. All model curves are normalized to 49 lenses, but their shapes 
are not tuned. 



